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. We study the conductivity of disordered zigzag graphene nanoribbons in the incoherent 

regime by using the Boltzmann equation approach. The band structure of zigzag nanoribbons 
contains two energy valleys, and each valley has an excess one-way channel. The crucial point 
is that the numbers of conducting channels for two propagating directions are imbalanced 
in each valley due to the presence of an excess one-way channel. It was pointed out that as 

, a consequence of this imbalance, a perfectly conducting channel is stabilized in the coherent 

regime if intcrvalley scattering is absent. We show that even in the incoherent regime, the 
conductivity is anomalously enhanced if intervalley scattering is very weak. Particularly, in 
the limit of no intervalley scattering, the dimensionless conductance approaches to unity 

, with increasing ribbon length as if there exists a perfectly conducting channel. We also show 

that anomalous valley polarization of electron density appears in the presence of an electric 
field. 

KEYWORDS: graphene nanoribbon, perfectly conducting channel, incoherent transport, 
J> ■ Boltzmann equation 
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CN | 1. Introduction 

(N : 

Graphene nanoribbons with zigzag edges (zigzag nanoribbons) have a unique band struc- 
| ture 1 ) which is not seen in ordinary quantum wires realized in semiconductor nanostructures. 

J> . Its band structure contains two valleys, called K_ valley and K + valley, well separated in 

•rH , 

^ ■ momentum space, and each valley has an excess one-way channel arising from a partially 

flat band. The crucial point is that the numbers of conducting channels for two propagating 
directions are imbalanced in each valley due to the presence of an excess one-way channel. 
This causes an unusual transport property of disordered zigzag nanoribbons in the coher- 
ent transport regime at low temperatures. Wakabayashi et al. 2 ' 3 ^ have pointed out that if 
impurity potentials are long-ranged and thus intervalley scattering is absent, one perfectly 
conducting channel without backward scattering is stabilized due to the imbalance between 
the numbers of conducting channels. This results in the absence of Anderson localization. 
Inspired by this observation, the statistical behavior of the conductance in disordered wires 
with channel-number imbalance has been studied extensively. 4 ~ n ) The presence of a perfectly 
conducting channel in such systems has been suggested by Barnes et al. 12 ' 13 ^ more than a 
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decade ago. However, since realistic systems has been lacking, detailed studies on this subject 
were not performed until recently. 

In addition to zigzag nanoribbons, a perfectly conducting channel appears in disordered 
carbon nanotubes. 14-16 ) In the case of carbon nanotubes, it is stabilized due to the two facts 
that the reflection matrix has the skew-symmetry and the number of conducting channels is 
odd. 16-18 ) Because the skew-symmetry plays a role only in the presence of the phase coherence 
of electrons, we expect that a perfectly conducting channel is fragile against dephasing. In- 
deed, Ando and Suzuura 16 ^ showed that a perfectly conducting channel and related anomalies 
disappear in the incoherent regime where the phase coherence of electrons is completely lost 
by strong dephasing. Recently, a quasi-perfectly conducting channel has been shown to appear 
in disordered graphene nanoribbons with armchair edges. 19 ) 

Note that the presence of a perfectly conducting channel in zigzag nanoribbons is con- 
cluded from the fact that the reflection matrix has a nonsquare form with dimensions 
(N + l)xJV due to the imbalance between the numbers of conducting channels, 2 ' 6 ' 12 ) where 
N + 1 (TV) is the number of incident (reflection) channels. From this fact we can show that 
one reflection eigenvalue becomes zero and thus one transmission eigenvalue becomes unity. 
This results in the presence of a perfectly conducting channel. Because the reflection matrix 
has no special symmetry in this case, we expect that a perfectly conducting channel is robust 
against weak dephasing. However, it should be noted that the reflection matrix itself is no 
longer well-defined in the incoherent regime because a single-particle description of quantum 
mechanical electron scattering is not justified under strong dephasing. That is, the above ar- 
gument cannot be applied to the incoherent regime, and therefore we cannot conclude the 
presence of a well-defined perfectly conducting channel. It is of interest and of significance to 
study whether anomalous electron transport properties arise in the incoherent regime. 

In this paper we show that even in the incoherent regime, the conductivity of disordered 
zigzag nanoribbons shows anomalous behaviors if impurity potentials are long-ranged and 
thus intervalley scattering is very weak. We introduce two parameters k and k' representing 
the strengths of intravalley scattering and intervalley scattering, respectively, and analytically 
obtain the conductivity as a function of k'/k by using the Boltzmann equation approach. 
Note that k' <C k if impurity potentials are long-ranged and k' « k if impurity potentials are 
short-ranged. We show that the conductivity is strongly enhanced if k'/k <C 1, while such an 
enhancement disappears in the case of k'/k = 1. Particularly, in the no intervalley scattering 
limit of k'/k — > 0, the dimensionless conductance approaches to unity with increasing ribbon 
length as if there exists a perfectly conducting channel. We also show that anomalous valley 
polarization of electron density appears in the presence of an electric field when k! <C k. We 
set h = 1 in the following. 
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Table I. The number N-r of right-moving channels and the number JVl of left-moving channels in 
each valley. 
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Fig. 1. The band structure of a zigzag nanoribbon consisting of 30 zigzag chains. The left and right 
valleys correspond to K_ and K + valleys, respectively. 

2. Analysis Based on the Boltzmann Equation 

We consider zigzag nanoribbons in which the K_ valley has N + m right-moving channels 
and N left-moving channels, and the K + valley has N right-moving channels and N + m 
left-moving channels (see Table I). Note that m = 1 in zigzag nanoribbons, while N can 
be controlled by carrier doping. The case of m > 2 may be realized in multilayer zigzag 
nanoribbons. 4 ^ Figure 1 shows the band structure of a zigzag nanoribbon consisting of 30 
zigzag chains. If the Fermi level is at the point indicated by the horizontal dotted line, the 
numbers of conducting channels are characterized by N = 2 and m = 1. Hereafter, we 
abbreviate the right-moving channels in the K± valley as and the left-moving channels 
in the K± valley as L^, and refer to R~ and L + as majority channels and L~ and R + as 
minority channels. Let e^ n ^. 

( £ Lnk) be tne energy of an electron state with wave number k in 
the nth right-moving (left-moving) channel of the K± valley. We introduce the corresponding 
distribution function g z Xnk -, where z = +, — and X = R, L. The group velocity v z Xnk is given 
by v z Xnk = de Xnk /dk. For simplicity we assume that v Rnk = vr and t>£ nfc = vl with v = vr = 

-VL > 0. 

We consider a long zigzag nanoribbon of length L placed along the x axis and apply a 
constant electric field E in the negative x direction. This field accelerates electrons in 
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while electrons in are decelerated. We express the distribution function as 

9xnk = fcv{e z Xnk ) - sign(v x )eE l z Xnk -^(e z Xnk ), (1) 

where /fd(£) is the Fermi-Dirac function and l Xnk , characterizing a deviation from the equi- 
librium distribution, is called the mean free path. This expression implicitly assumes that 
the applied electric field increases (decreases) the population in (L 1 * 1 ). The distribution 
function obeys the Boltzmann equation, 20 ) 

~df~ ~da~ ~dk~ ~ ^ W Xnk,X'n'k' [9 X'n'k' ~ 9Xnk ) > ( 2 ) 

X'n'k'z' 

where W z Xnk x , n , k , is the scattering probability between the state with {X'n'k'z'} and that 
with {Xnkz}. The scattering probability is expressed as 



Xnk,X'n'k' 

with 



W%* u v ,„,u, = ^M z / nk X , nlk ,8 (e z x , n , k , - e z Xnk ^j (3) 



M 2 ' 2 

lvl Xnk, X'n'k' 



u Xnk,X'n'k' 



(4) 




where U x z nk x , n , k , is the matrix element of impurity potential and (• • • ) indicates the ensemble 
average over impurity configurations. For simplicity we assume that the scattering probability 
is determined by only the valley indexes z and z', and does not depend on details of initial 
and final states. Hence M x '^ k x , n , k , = M and M x ^ k x , n , k , = M', where M and M' describe 
intravalley scattering and intervalley scattering, respectively. The Boltzmann conductivity a 
is obtained as 

a = 1 Z> \ Vx \ l Xnk I Qp( £ Xnk, 

N N N+m \ 

Rn + E / Ln + E / «n+ E l tn) > ^ 
n=l n=l n=l / 

where l Xn is the mean free path at the Fermi level. Here and hereafter the spin degeneracy is 
ignored. 

To uncover anomalous features of zigzag nanoribbons in an electric field, we focus on 
steady states with spatial uniformity. According to the assumptions stated above, the mean 
free path becomes independent of n. We thus set l z Xn = l x , where k is also dropped since we 
are interested in electron states near the Fermi level. Furthermore we can set I = = l£ 
and I' = Zl = from the symmetry of the band structure, where Z and I' correspond to the 
majority and minority channels, respectively. This ensures charge neutrality of our system. 
Substituting eq. (1) into eq. (2) and ignoring the t- and x-dependent terms, we obtain two 
equations for Z and I'. In terms of the parameters k and k' defined by 

LM (to 
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(7) 

the resulting equations are expressed as 

1 = [Nk + (SN + 2m)K'] I + N(k - k')1', (8) 

1 = [(iV + m)« + (3iV + m)«'] I' + (AT + m)(/c - (9) 

Note that k and k' represent the strength of intravalley scattering and that of intervalley 

scattering, respectively. We easily obtain 

I f AN m\ 
1 = ^T-ZZv>. -7—7, + 77 > (10) 



2(2iV + m) 2 \k 

= 1 / 4(iV + m) _ m\ 

2(2N + m) 2 \ k + k' k'J ' 1 j 

Substituting eqs. (10) and (11) into eq. (5), we obtain the Boltzmann conductivity as 

e 2 1 f8N(N + m) m 2 ' 



2vr (2N + to) 2 \ n + n' + k' ) ' ^ 
These results show that anomalous features arise in the case of to / 0. Firstly, we observe from 
eq. (11) that l£ and ^ become negative when k' is sufficiently small. That is, g^ nk increases 
and g^ nk decreases. This indicates that the distribution functions for the minority channels 
shift to the direction opposite to that assumed in eq. (1). It should be emphasized that because 
9xnk i ncreases an d 9xnk decreases regardless of the propagating directions, electron density 
is polarized between the two valleys in the presence of an electric field. Secondly, we observe 
from eqs. (10) and (11) that all the mean free paths diverge in the limit of k' — > as in 
ballistic quantum wires. Accordingly, the Boltzmann conductivity also diverges. 

To gain an insight into these anomalous features, we consider the Boltzmann equation 
including the t-dependent term in the limit of k' — > 0. The two valleys are completely decoupled 
in this limit, so we treat only the AT_ valley. The Boltzmann equation yields 



ldl 



V 



R 

dt 



Nk{I r +17), (13) 



1 dl~ 

~if = ^-(N + m)K(l R + l7). (14) 

In the right-hand side of each equation, the first term represents the acceleration induced by 
an electric field, while the second term describes the deceleration due to backward scattering. 
These equations indicate that the electric field increases l R and 17] in the same rate but the 
deceleration effect for l R is by factor N/(N + m) smaller than that for /£" due to the imbalance 
between the numbers of conducting channels. Therefore, the acceleration and the deceleration 
are never balanced in the case of m ^ 0, and l R and 17] do not converge. The acceleration 
dominates the deceleration in R~ (i.e., majority channels) while their relative strengths are 
interchanged in L~ (i.e., minority channels). This indicates that l R — > +oo and ~~ ^ — 00 • 
Similarly, we observe that l£ — > +oo and l R — > -co. We understand that the anomalous 
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increase of / (i.e., / — > +00) is caused by excess electric field effect while the anomalous 
decrease of I' (i.e., I' — > — 00) is caused by excess backward scattering. From the argument 
given above, the characteristic behavior of electrons in the limit of k! — > with m / is 
summarized as follows. Firstly, I — > +00 implies that the behavior of electrons in the majority 
channels is similar to that in ballistic quantum wires in spite of the presence of backward 
scattering due to impurities. Secondly, V — > — 00 implies that the behavior of electrons in the 
minority channels is also similar to that in ballistic quantum wires but the electric field effect 
arises in the direction opposite to that naturally expected. That is, electrons behave as if they 
are free from impurity scattering but an electric field affects electrons in the minority channels 
as if its sign is reversed. 

The divergence of the mean free paths must be removed if we correctly take account of 
the fact that electrons are accelerated by an electric field only in a finite region. Therefore we 
reconsider the Boltzmann equation including the x-dependent term. Suppose that a constant 
electric field E is applied only in the region of —L/2 < x < L/2. We consider that incident 
electrons from left (right) are described by the equilibrium distribution at x = —L/2 (x = 
L/2). This indicates that 21) 

/ R (-L/2) = Z+(-L/2)=0, (15) 
£(L/2) = £(L/2) = 0. (16) 
From the symmetry of the band structure, we assume that 

l(x) = = l+(-x), (17) 
l'(x) = l+(x) = lZ(-x). (18) 
The Boltzmann equation for steady states yields 

+ 1 = [ Nk+ ( 2N + m ) K '] l ( x ) + Nk1 '(~ x ) 

+ (N + m)n'l{-x) - Nk'I'(x), (19) 

dl'(x) 

dx +1= + + (y2N + + ( N + m W(~ x ) 

+ Nk'I'(-x) - {N + m)Kl(x). (20) 
We solve these coupled equations under the boundary condition, eqs. (15) and (16), and obtain 



L L 
l{x) = x + - + - 



l'{x) = x + - + - 



m(n - k) (c(x) - c ) H = (d(x) - do) 



- 2(2iV + m) 2 (n + k')k'cq(x + ^) 
m(/c - k') (c(x) - c ) + 171 ( K K ) (d(x) - d ) 

yJOL 



(21) 
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(22) 



- 2(2iV + m) 2 (K + k')k'cq(x + ~; 



where a = (k + K,')[m 2 n + (8iV 2 + 8Nm + m 2 )*!] and 

. , (2N + m)(K + k') , . _ . , . _ . 

c(x) = -= cosn(va^) — smn(vax), (23) 

ya 

.. , (2N + m)(K + k!) . . . r- , . . ,- . 

d(a;) = j= sinh(Vax) + cosh(V«^)- (24) 

The constants in the above equations are given as follows: cq = c(— L/2), do = d(—L/2), and 

f\ \ k 2 - k' 2 

S = 4(2iV + m)/e' 1 + -(2JV + m)(/e + k')L c + 2m 2 ^d . (25) 

V 2 J y/a 

We now obtain the Boltzmann conductivity cr(x) and the two-terminal conductance G. 

Equation (5) yields 

e 2 

<*(x) = 7T- [(^ + m ) (K^) + K-^)) + ^ (*'(*) + «'(-a:))] • (26) 

Substituting eqs. (21) and (22) into eq. (26), we find that a(x) is independent of x and is 
given by 

= e 2 L 2ac 

2vr (k + k')S' 1 ^ 

This ensures current continuity in our system. The conductance G = a/L is expressed as 

e 2 [m 2 n + (8A^ 2 + 8Nm + m 2 )/?'] c 

G = — 2 — • (28) 

2vr 2 (2iV + m)K' [l + ±(2iV + m){n + k')L] c + m 2 ^^d 

Let us consider the no intervalley scattering limit of k' — > 0, at which we obtain 

e 2 / 2Nm \ 

2^\ m+ (N + m)e™ L - N J ' ( } 

We observe that the dimensionless conductance g = (2n/e 2 )G behaves as g — > m in the 
limit of kL —t- oo as if there exist m perfectly conducting channels. This indicates that the 
conductance in the limit of k' — > is drastically enhanced and greatly deviates from Ohm's 
law. We interpret it on the basis of the argument given below eqs. (13) and (14). We show that 
very anomalous electron distributions are induced by an electric field in the limit of k' — > 
with m 0. These distributions imply that electrons behave as if they are free from impurity 
scattering as in ballistic quantum wires but the electric field effect on L~ and R + arises as if 
the sign of an electric field is reversed. This observation indicates that the contribution to g 
from R~ and L + is N + m while that from L~ and R + is negative and is given by —N. By 
summing these two contributions we obtain g = m. This roughly accounts for our result in 
the limit of kL — > oo. 

Before considering the case of k' > 0, we briefly examine the two limiting cases in which 
Ohmic behavior of G is expected. In the absence of one-way excess channels (i.e., m = 0), 
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Fig. 2. The dimensionless conductance as a function of kL for k' /k = 0.0001, 0.01, and 1.0 in the 
case of N = 2 and m = 1. 



eq. (28) is reduced to 

G " 2vrl + iV(re + re')L' {6U) 
which obeys Ohm's law when iV(re+re')L 3> 1. We next consider the case of re' = re in which the 
two valleys are completely connected and the imbalance between the numbers of conducting 
channels disappears. In this case, eq. (28) is reduced to 

e 2 2N + m 

2vrl + {2N + m)nV { ' 

which also obeys Ohm's law when (2N + m)reL 3> 1. 

We now consider how the behavior of the conductance is affected by re' in the case of 
m ^ 0. We numerically obtain the dimensionless conductance g = {2it/e 2 )G as a function 
of kL or re'/re. Figure 2 shows g as a function of kL for re'/re = 0.0001, 0.01, and 1.0 in the 
case of N = 2 and m = 1. We observe that 5 is nearly unity even when kL = 20 in the case 
of re'/re = 0.0001. Figure 3 shows g as a function of re'/re for reL = 10 and 20 in the case of 
N = 2 and m = 1. Since the behavior of 5 for small re' is of interest, we treat only the region 
of re'/re < 0.01. The result for the case of N = 3 and m = is also displayed for comparison. 
We observe that in the case of N = 2 and m = 1, g approaches to unity with decreasing re'/re, 
while g in the case of N = 3 and m = is nearly independent of re'/re. Furthermore, 5 in the 
case of N = 2 and m = 1 is much larger than that in the case of iV = 3 and m = in spite 
of the fact that the total number of conducting channels in the latter case is by one greater 
than that in the former case. This indicates that the conductance is strongly enhanced in the 
case of m ^ if re' is sufficiently small and thus the intervalley scattering is very weak. 

3. Summary and Discussion 

In summary, we have studied the conductivity of disordered zigzag graphene nanoribbons 
in the incoherent regime by using the Boltzmann equation approach. We have shown that 
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Fig. 3. The dimensionless conductance as a function of k'/k for kL = 10 and 20 in the case of TV = 2 
and to = 1, and that in the case of N = 3 and m = 0. 

the Boltzmann conductivity of zigzag nanoribbons is anomalously enhanced if intervalley 
scattering is very weak. Particularly, in the limit of no intervalley scattering, the dimensionless 
conductance approaches to unity with increasing ribbon length as if there exists a perfectly 
conducting channel. We have also shown that anomalous valley polarization of electron density 
appears in the presence of an electric field if intervalley scattering is very weak. 

We here consider how the dimensionless conductance g of disordered zigzag graphene 
nanoribbons behaves with increasing temperature, focusing on the effect of dephasing. For 
simplicity we restrict our consideration to the most interesting limit of no intervalley scat- 
tering. At zero temperature where the effect of dephasing can be ignored, we observe that 
one channel becomes perfectly conducting and the other channels are subjected to Anderson 
localization. 2 ' 3 ) Hence g decreases exponentially toward unity with increasing ribbon length. 
The effect of dephasing becomes stronger with increasing temperature. However, as long as 
dephasing is weak, we expect that one perfectly conducting channel survives as suggested in 
§1 and therefore g behaves as in the zero temperature limit. Our result eq. (29) indicates 
that the behavior of g is almost unchanged even if temperature is further increased and the 
dephasing effect becomes strong. This should be contrasted to the behavior of g in disordered 
carbon nanotubes in which a perfectly conducting channel and related anomalies are strongly 
suppressed by dephasing. 16 ) We conclude that the temperature range in which the anomalous 
enhancement of the conductivity can be observed is much wider in zigzag nanoribbons than 
in carbon nanotubes. 
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